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The equation of state for two flavor QCD at finite density * 
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We discuss the equation of state for QCD at non-zero temperature and density. We present results of a 
simulation for QCD with 2 flavors of p4-improved staggered fermions. Derivatives of In Z with respect to quark 
chemical potential Hq up to fourth order are calculated, enabling estimates of the pressure, quark number density 
and associated susceptibilities as functions of Hq via a Taylor series expansion. We also discuss the radius of 
convergence of the expansion as a function of temperature. It is found that the fluctuations in the quark number 
density increase in the vicinity of the phase transition temperature and the susceptibilities start to develop a 
pronounced peak as Hq is increased. This suggests the presence of a critical endpoint in the (T, /Xq) plane. 



1. Introduction 

In the last few years, remarkable progress has 
been made in the theoretical study of QCD ther- 
modynamics with small but non-zero quark chem- 
ical potential ^iq. In particular, numerical studies 
of the phase structure in the {T,fiq) plane [1-3] 
provide important information for heavy-ion col- 
lision experiments. One of the most interesting 
investigations is finding the critical endpoint, at 
which the crossover phase transition line becomes 
first order. It might be possible to detect the end- 
point experimentally via event-by-event fluctua- 
tions, since the number density fluctuation should 
be large around the endpoint [4] . The fluctuation 
should be in proportion to the susceptibility of 
the quark number density, hence the study of the 
equation state at non-zero baryon number density 
is very important [5,6]. 

The quark number susceptibility is given by 
the second derivative of pressure with respect to 
Several studies of the quark number suscep- 
tibility have been performed at Hq — [7], and 
measurements of pressure and energy density at 
7^ were done by [9] using the reweighting 
method, which allow us to investigate thermo- 
dynamic properties at fiq ^ 0. This approach. 
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however, does not work for large fiq and large 
lattice size due to the sign problem. Here, we 
adopt the following strategy. We compute the 
derivatives of thermodynamic quantities such as 
pressure, quark number susceptibility etc., with 
respect to fiq a.t ^q — 0, and determine the Tay- 
lor expansion coefficients in fiq [2,8], in which the 
sign problem does not arise. 

We discussed last year the relation between 
the line of constant pressure (energy density) and 
the phase transition line calculating the second 
derivative of pressure [2]. In this study [10], we 
proceed to fourth order derivatives. Then, we can 
for the first time investigate the /ig-dependence 
of the quark number susceptibility near fiq ~ 0. 
(Because the pressure is an even function of fiq, 
the fiq-teim is leading and /i^-term is the next 
to leading, and in fact only these two terms are 
non-zero in the high temperature limit.) By esti- 
mating the change of the susceptibility, we discuss 
the possibility of the existence of the critical end 
point in the phase diagram of T and fiq . 

2. Taylor expansion in fiq 

Pressure is given in terms of the grand partition 
function Z(T,V,fiq) by p/T^ = {l/VT^)liiZ. 
However, the direct calculation of \nZ is diffi- 
cult, hence most of the work done at /x^ = for 
the calculation of pressure is done by using the 
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Figure 1 . CocfRcients of Taylor expansion, C2 (up- 
per) and C4 (lower). Tq is at /Zg = 0. 

integral method [11,12], where the first deriva- 
tive of pressure is computed by simulations, and 
the pressure is obtained by integration along a 
suitable integral path. For fiq 0, direct Monte 
Carlo simulation is not applicable; in this case we 
proceed by computing higher order derivatives of 
pressure with respect to fiq/T a,t iJq — 0, and then 
estimate p(/Ug) using a Taylor expansion, 
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where c„ = (l/n!)5"(p/T4)/a(/x,/T)"|^,=o. 
These derivatives can be compiitcd by the ran- 
dom noise method, which saves CPU time, and 
also it can be proved that the odd terms are ex- 
actly zero. Furthermore, wc do not need simula- 
tions at (T, Hq) = (0, 0) for subtraction term of p 
to normalize the value of p, since the derivatives 
of p|T=o,Mg=o with respect to ^q are, of course, 
zero. This also reduces CPU time. 

The quark number density Uq and quark num- 
ber susceptibility Xq are obtained by the deriva- 



tives of pressure; 

nq _ d{p/T^) 
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We compute the pressure up to /Xg-term using 2 
flavors of p4-improved staggered fermions [13] at 
ma = 0.1 on a 16'^ x 4 lattice. Then the quark 
number density and quark number susceptibility 
are obtained up to 0{p?q) and 0{p.'^), respectively. 
The details of the simulations arc given in [10]. 
In Fig. 1, we plot the data for C2 (upper) and 
C4 (lower). Both of them are very small at low 
temperature and approach the Stofan-Boltzmann 
(SB) limit in the high temperature limit, as we 
expected. The remarkable point is a strong peak 
of C4 around Tc- 

We can understand this peak by the follow- 
ing two issues. One is a prediction from the 
hadron resonance gas model [14], which is an ef- 
fective model of the free hadron gas in the low 
temperature phase. The model study predicts 
Ci/c2 = 0.75 and our results arc consistent with 
this prediction for T < Tc, in fact, as T increases 
C4/C2 remains constant until T ^ whereupon 
it approaches the SB limit. 

The other point is from a discussion of the 
convergence radius of the Taylor expansion. We 
expect that the crossover transition at ^iq = 

changes to first order transition at a point 
Hq/Tc ~ 0(1) [1,15]. Then, the analysis by 
the Taylor expansion must break down in that 
regime, i.e. the convergence radius may be 
smaller than the critical endpoint. We de- 
fine estimates for the convergence radius by 
Pn = \/\cn/cn+2\- We computc po and p2 from 
p/T^{Qi), C2 and C4. It is found that both pQ and 
P2 arc quite large at high temperature as expected 
from SB limit, ~ 2.01, ~ 4.44. On the 
other hand, around these are 0(1), since C2 
and C4 arc of the same order, so that our results 
around Tc suggest a singular point in the neigh- 
borhood of fJiq/Tc = 1. 

3. Equation of state at pq ^ 

Next, we calculate pressure and quark num- 
ber susceptibility in a range of < PqjT < 

1 which is within the radius of convergence 
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Figure 2. Difference of pressure from = as a 
function of T for each fixed [iqjT . 




Figure 3. Quark number susceptibility as a func- 
tion of T for each fixed [iqjT . 

discussed above, using the data of C2 and 
C4; A(p/T4) ^ v{T,ii,)lT^ - v{T,^)lT^ = 
C2{i^jTf + Ci(pjTY + 0(/i6), and xjT^ = 

2c2 + 12c4(m,/T)2 + 0(;,4)^ 

We draw A{p/T^) for each fixed fig/T in Fig. 2 
and find that the difference from p|^^=o is very 
small in the interesting regime for heavy-ion colli- 
sions, fig/T « 0.1 (RHIC) and Hg/T « 0.5 (SPS), 
in comparison with the value at fXq = 0, e.g. SB 
value for 2 flavor QCD at fig = 0: p^^ /T^ ~ 4.06. 
The finite density effect of pressure at fiq/T = 0.1 
is only 1%. This conclusion is consistent with that 
inferred from just the 0{fi^) term [2]. Also, the 
result is qualitatively consistent with that of [9] 



obtained by the reweighting method. 

Figure 3 is the result of Xq/T"^ for fixed fiq/T. 
We find a pronounced peak for mq/T > 0.5, 
whereas Xq does not have a peak for fiq = 0. This 
suggests the presence of a critical endpoint in the 
{T,fiq) plane. 

This discussion can be easily extended to the 
charge fluctuation xc- The spike of xc at Tc is 
weaker than that of Xq, which means the increase 
of the charge fluctuation is smaller than that of 
the number fluctuation as fMq increases [10]. 
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